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Abstract 

In this paper we investigate the Q-ball Ansatz in the baby Skyrme 
model. First, the appearance of peakons, i.e. solutions with extremely 
large absolute values of the second derivative at maxima, is analyzed. It 
is argued that such solutions are intrinsic to the baby Skyrme model and 
do not depend on the detailed form of a potential used in calculations. 
Next, we concentrate on compact non spinning Q-balls. We show the 
failure of a small parameter expansion in this case. Finally, we explore 
the existence and parameter dependence of Q-ball solutions. 

1 Introduction 

The baby Skyrme model is an important member of the family of nonlinear 
field theories. Although it was proposed as a simplified version of the Skyrme 
model, it is used to model some real- world systems, see [T], [5] and references 
therein. From the formal point of view, the theory describes dynamical maps 
between two spheres S2. This gives rise to a non-trivial topological structure 
and most efforts were done in this domain. However, in this paper we shall 
put the topological questions aside since the trivial topological sector is also 
quite interesting, see [3J. If a customary field potential is chosen, there is an 
additional, apart from the topological charge, quantity untouched by the time 
evolution. In the scalar field formulation it is simply the U(l) charge. Its pres- 
ence validates the question about existence of Q-balls, solutions known in many 
theories to minimize the energy for a given value of charge. 
The present work aims at clarifying some findings described in [TJ. In that 
paper C. Adam et al. report on existence of both spinning and non-spinning 
compact Q-balls in the baby Skyrme model with a V-shaped potential. Some 
of these solutions, called peakons, have a peculiar shape - the value of the sec- 
ond derivative at their maxima becomes extremely large. These observations 
come from numerics and some of them lack satisfactory analytical explanation. 
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In the present paper we shed some light on them. First, we identify a struc- 
ture in equation for the Q-ball profile function responsible for the appearance 
of peakons. Second, we describe qualitatively the existence domain of the non- 
spinning Q-balls in the theory. On the way, difficulties with the customary small 
parameter expansion in theories with non smooth field force are discussed. 
We consider a sharp potential associated with the signum-Gordon model (see 
e -g- ID) giving rise to compact solutions and requiring a matching procedure. 
Some objections to the potential has been risen in [5]. On the other hand, 
the article [5] may be read as an apologia for such potentials. The solutions 
described in the present paper are legitimate weak solutions to a differential 
equation, a discussion of this fact may be found in [7]. Another attitude to- 
ward the potential is also possible. Then it holds as an idealization, where 
the close-to-the- vacuum effects are neglected. The results in [5] suggest that a 
regularization of the potential would not change seriously the below described 
solutions. 

The paper is organized as follows. After a short introductory section we dis- 
cuss the well-posedness problem of the equation for the Q-ball profile function. 
This issue is intimately connected to the appearance of peakons. In the next 
section we give an account of our futile attempt to formulate a small parameter 
expansion for Q-balls in the theory at hand. A natural and well working ap- 
proximation is introduced in the section 5. It makes it possible to characterize 
the Q-ball dependence on parameters present in the equation. The last section 
contains our conclusions and remarks. 



2 General settings 

Let us consider the following Lagrange density function: 



(d^d^®) 2 - d^d^Qd^d"® |$| 

: (i + |$| 2 ) 2 _8/? (f + |i>| 2 ) 4 A V /TT^| 



where $ is a complex scalar field, $ its complex conjugation and |<I>| 2 = $$; 
j3 > and A > 0. The above form of the Lagrangian in three (2 + 1) dimensional 
Minkowski space-time defines the baby Skyrme model in the scalar field formu- 
lation, for details see pQ. The last term in the Lagrangian is a field potential. 
In this paper we use the V-shaped one following [i,. The charge Q coming from 
the U(l) phase invariance is given by the formula 



Q = Im Id x 



$9 n $ 16,3$ , - - - -s 



(i + \m 2 (i + i*i 2 ) 4 

"(2) 

Let us stress that this quantity is not related to the topological structure of the 
theory 

We shall consider solutions of the model set by ([1]) found with help of the 
spinning Q-ball Ansatz: 

$(x) = e i(a " + "^/(r), 

where to > 0, (r, </>) are radial and angle coordinates, / is a real valued function 
and finally n is an integer. In [T] it is pointed out, that for small amplitudes 
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the theory reduces to the signum- Gordon model. This fact suggests writing the 
equation in a form corresponding to the on-shell scaling symmetry present in the 
signum-Gordon theory. Rescaled variable y — ujr and function Xg(y) — 8uj 2 f(r) 
should be appropriate and lead to some simplifications. In this setting we obtain 
the equation 
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+ signig)^! + K 2 g 2 , (3) 



where k — \/8u> 2 and 7 = A/3. Thus, the equation depends on two combinations 
{7, k} of the original three parameters {A, /3, uj}. The function sign(g) = 1 
if g > and sign(g) = if g = 0. We shall be interested only in nontrivial 
solutions with g > 0, hence we omit the function sign(-) in what follows. The 
solutions of physical importance satisfy the boundary conditions ensuring finitc- 
ness of the energy and the charge. So, before we proceed with the analysis of 
the equation ((3]), let us write down the expressions for the energy and the U(l) 
charge in terms of the rescaled variable y and function g 



E = 2ttk z 



dy y 
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(5) 



In what follows we vary k while 7 is kept constant. This choice of the control 
parameter is rather natural. For a given model the parameters A and ft (thus 7) 
are set. The parameter k depends also on ui that characterizes initial data. 
Hence, it may be any real number. 



3 Appearance of peakons 

Henceforth we consider the equation ([3]) for n = 0. We rewrite it in a form 
convenient to further analysis 

1 K 9 2 \ n , 1 (« IKg 2 \ , 

' 7. 1 - M , )9 + (6) 



(1 + K 2 g 2 ) 2 / y \ (1 + n 2 g 

9 (i + .v )3 - 9 + 1 1 "tt^vJ 3 = Vi + kV - 

The left-hand side of the above equation does not depend on the specific po- 
tential and is peculiar to the Q-ball Ansatz in the baby Skyrme model. It is 
a differential equation of the second order, hence for given yo > 0, g(yo) and 
g'(yo) we should be able to set uniquely g"(yo)- But this is not true for some 
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values of g. The above form of the equation makes this evident: the factors 
multiplying g" and g' may be simultaneously set to zero. We shall argue below 
that this lack of well-posedness gives rise to peakons. To this end we examine 
solutions of (O in a vicinity of g r , where g r satisfies 

1 - ^ = 0. (7) 

As for the number of roots of the above relation and their dependence on the 
parameters k and 7 see figure [TJ For convenience we introduce a new function 
£ {y) = .9(2/) — 9r- In the vicinity of e = the following equation holds 

ee" + -ee' + a l {e'f = a 2 , (8) 
V 

where a's are constants depending on g r , 7 and k. The equation is obtained 
from ((5]) by expansion around g r of the coefficients multiplying g", g' , g' 2 and 
the terms depending on g only. The leading terms are retained. The value of 
ai is entirely determined by the baby Skyrme model. It is equal to 1/2 for 
< k < K m ax, for n m ax it vanishes. In what follows we shall not put the actual 
value of a\ to make our considerations as general as possible. We shall see that 
the exact value of the parameter does not play an important role. The potential 
partially determines the value of a 2 - The equation is valid unless one of 
its parameters {ol\, a 2 ) vanishes or explodes. For a given 7 it happens for few 
discrete values of k. Thus, the equation ([8]) is a generic one and we concentrate 
on it. 

Let us state two immediate observations about the equation (|S]). First, the 
reflection e —> — e does not alter the equation. Second, rescaling of the variable 
V \f\& 2 \\) makes the constant term on the right-hand side in ([5]) equal to ±1. 
Hence, a 2 sets the scale but its exact value is not crucial. The value of oi\ 
cannot be altered by rescaling of the function s or variable y. 
Let us now analyze the case e' = 0. Then we have e" — a 2 / £• It means, e is 
repelled from the value e = if c*2 > and attracted to it if 02 < 0. To be 
more concrete, consider e < 0. If 02 > then e has a maximum, otherwise 
minimum. Note, that the smaller absolute value at the extremum, the sharper 
the extremum is. This remark is valid even if y = is considered. 
Let us note, that the equation ([S]) has a simple solution regular for y — 0. It 
reads 

e(y) = ±\ 77-7— ;y- 

V (1 + 011) 

The equation (JSj) is not analytically solvable in general. A great deal of infor- 
mation about the solutions may be inferred from the equation. To this end the 
substitution 

e(y) = \w(y)\^ (9) 

is very useful. The function w has only nonnegative values. The equation ([8]) 
turns into 

1 , d ( (1 + tti) 2 ^\ nn , 

w + -w = -— -a 2 — ■ (10) 

y aw \ 2«i J 

The above form of the equation makes its mechanical interpretation natural. 
The equation (|10[) is seen as a Newton equation for a particle in a potential sub- 
ject to a friction (the term with the first derivative). Assume 2a\/(ai + 1) > 0. 
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Figure 1: The structure of the roots of the relation ([7]) for given 7. The maximal 
value of k for which a real solution exists is K ma x = 7/4. A minimal value g r for 
a given 7 is equal to 16/(\/277) and is reached for k = 37/I6. For the value 
Ki it holds OI2 = 0; the blue line (on the right from Ki) marks 012 > and the 
for red line (on the left from ki) ot2 < 0. For small 7 the two colors meet on 
the upper branch, for large values of 7 on the lower one. 



Then, it is the sign of the coefficient a^/ai that decides about the qualitative 
behavior of the solutions. If it is negative, the particle is trapped in a po- 
tential well. Any solution eventually reaches the zero value then. This value 
is approached with a nonvanishing first derivative. Hence, the solution of ([8]) 
approaches zero as e ~ 

{yi - yY^ ai+1 \ This gives rise to the exploding first 
derivative of e when y —> y±. Thus, there is no unambiguous way to continue the 
solution for y > y±. Let us briefly interpret this finding in terms of the original 
equation ©. In a neighborhood of g r any solution of (JBJ coincides with the 
corresponding solution of @. Thus, any solution coming close enough to g T has 
a dead-end since it eventually reaches the value g r with infinite first derivative. 
Now we discuss the case of positive otz/otx in equation ([10]) . Then, the po- 
tential has a global maximum for w = 0. A fictitious particle can climb the 
potential maximum if it is given enough (mechanical) energy. This gives rise to 
the exploding first derivative of e as described above. However, this is not the 
only possible scenario. The particle may not have enough energy to reach the 
potential maximum. Then, it climbs for some time (y is seen in this context 
as time), reaches its maximal accessible potential energy and finally dips. Of 
course, up to the initial data, the particle may just decrease its potential energy 
without approaching the maximum. A corresponding solution reads asymptot- 
ically w ~ y ai , hence e(y) ~ y. 

There is also a solution separating the two just described families of solutions. 
It is a solution grazing the value w = at a point y±. The function w behaves 
then like (yi — y) ai+l , what gives e(y) w ±y/a2/ati\yi —y\. Precisely, the equa- 
tion © is ambiguous for e = and the behavior e ~ ±(y — yi) is also correct. 



■5 



9 




9r 



6.4 



6.3 



6.2 



6.1 



0.32 



0.37 0.42 0.47 y 



0.415 



0.421 



0.427 



0.433 



y 



Figure 2: Numerical solutions of the equation ([5]) with 7 = 0.5 and K = 0.07. 
This implies g r — 6.427 and U2 — 15.34. Initial data are set for y = 0.4 and 
g = 6.2. On the right panel the derivative g' of the solution approaching closer 
to g r is 5.745, the other solution starts with g' = 1.5. The solution on the right 
panel is solved with g' = 6. 



Let us turn to the equation ©. Again, we concentrate on a strip around g r , 
where the equation ^ turns effectively into ©. Consider solutions entering 
this strip (thus the value of g is set) at a given point. The further evolution 
depends on the derivative of g. If it is small, the function has an extremum and 
leaves the strip. The bigger absolute value of the derivative, the closer to g r the 
function g approaches. For a critical value of the derivative the grazing solution 
appears. For larger values, the solutions touch g r with infinite first derivative. 
Figure [5] illustrates the above discussion graphically. 

In light of the above analysis the appearance of peakons is just a question of 
amplitude. If a Q-ball profile function comes close to a root of the relation (|7|) , 
then a peaklike solution is expected (if aij(X\ > 0). Let us stress once more 
that the above described structure does not emerge due to the V-shaped poten- 
tial. It is inherent in the Q-ball Ansatz in the baby Skyrme model. It would 
be interesting to see peakons in a broader perspective: what is their role in the 
whole model and what kind of nonlinearity makes them possible. 
The peakons have been reported in the case of the spinning Q-balls, see pQ. In 
this case the equation for the profile function is a bit more complicated. We 
have decided to work out the problem with n — as it seems both generic and 
technically simple. What is more, some peakonlike nonspinning Q-balls are also 
present in the model. We shall see this in the fifth section. In the next section 
some technical aspects peculiar to V-shaped potentials are described. 
Finally, one more remark is in order. An analysis of equation ^ for large values 
of the function g shows, that solutions explode to infinity for finite y. We are 
not about to delve into this here, though. 
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4 The case k = and power expansion 



Now, we are ready to look at Q-balls. We begin the analysis of the equation © 
with investigation of the case k = 0. The equation has then the following form 

g" + -g' + g = sign(g), (11) 
y 

coinciding with the equation for the Q-ball profile function in the signum-Gordon 
model. The relevant solutions in this model are well-known, see [4] and [8] . One 
can expect that the Q-ball solutions of ([5]) may be expressed in the (traditional) 
power series of the parameter k: g = go + ngi + K 2 g 2 + ■ • • for k ranging from 
zero to some maximal value. Note, that we can consider initial data with uj as 
large as needed. This means, n may be arbitrarily close to zero. The boundary 
conditions supplementing both the equations (j6]) and (fTTj) are as follows. To 
ensure regularity of the solution at the origin it is necessary that g'(0) = 0. The 
energy and charge are finite if g{y) — > if y — > oo. The nonsmooth field force 
makes the matching procedure indispensable, i.e. we need to find a point yo > 0, 
for which g(yo) = g'(yo) = 0. For y < yo the function g is strictly positive, for 
y > yo the profile function vanishes. Let us now try to find a solution of the 
equation ([6]) in the following form 

g = .9o + + K 2 g 2 + ..., 

where go is a solution of (fTTj) . To this end the above Ansatz is plugged into the 
equation. In the first order we obtain the relation 

," , 1/ , ( 9o9o . /2 2 //2 



<?i + ~g\ + .91 = -7 I — — + 5oSo + SoSo ) • ( 12 ) 



It is solved using the Green function technique: 



gi{y)= dssG(y,s)F(g (s))+A 1 J (y)+A 2 Y (y), (13) 
Jo 

where F denotes the r.h.s. in the equation (|12j) . Jo an d Yq are Bessel functions 
of the first and second kind solving the homogeneous part of the equation (TT2"j) . 
Note, that J o (0) > and J o (0) = 0, but Y (y) -> oo when y -> 0. Thus, due 
to the boundary conditions A 2 — and A\ is a free parameter. We use the 
following Green function 

n , Ms)Yo(y)~Y ( S )Jo(y) 
[V > ) y{Y>[y)Jo{y)-J>MY {y)y 

The integration in fj 13[) gives a function behaving like y 2 In y at the origin. This 
Green function does not assume any condition at any other point y > 0. This is 
why the freedom in adding a term proportional to Jo is left (for further details 
consult [5]). This freedom is necessary to make the function go + ngi satisfy the 
boundary conditions. One can choose a different Green function with explicit 
dependence on the (unknown) matching point, see [1]. However, it would not 
be helpful. To proceed, it has to be specified how to deal with the boundary 
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conditions. We require that in any order of the power expansion the following 
conditions be satisfied: for any m there be a point y m such, that 



go{ym) + Kgi{Vm) + • • ■ + K m g m (y m ) = 0, 
g'oiVm) + Kg'iiVm) + ■■■ + K m g' m {y m ) = 0. 



To this end the coefficient A\ (and analogously in higher powers) has to vary 
with k, making 171 K-dependent. Thus, the boundary conditions make the ex- 
pansion fail. Another ambiguous decision to be made is how to understand go 
for arguments larger than ?/o- At this point the equation with nonsmooth force 
term is not uniquely determined, see [8]. To sum up, the power series expansion 
is just an abuse of notation. Fortunately, we can resort to a more successful 
approximation. 



In the case of an analytical potential giving rise to Q-balls there is an intuitive 
tool to cast a light on the relevant equation: the mechanical analogy. As de- 
scribed in the third section the equation for the profile function is then viewed 
as a Newton equation for a particle moving in an effective potential and subject 
to a friction, see e.g. [9]. A Q-ball profile function is then usually a solution 
interpolating between three static solutions: one of them is stable, two are not. 
This picture is simplified in the signum-Gordon model. If sign(g) in the equa- 
tion ([lljl is taken to be +1, then the resulting equation has a simple structure: 
there is only one constant and stable solution. Any other solution may be viewed 
as its linear perturbation. It is the role of the boundary conditions to select a 
perturbation interpreted as a Q-ball solution. We shall adapt this idea to the 
baby Skyrme theory. Obviously, equation is a nonlinear one and any lin- 
earization has a limited range of validity. 

First let us find constant solutions of equation ©. They satisfy the following 
relation 



It turns out, it has zero, one, two or three real roots. Their structure illustrates 
figure [3] There are three branches of solutions: c\, C2 and C3. In what follows 
Ci denotes a branch of roots in general or a specific root from the branch Cj. 
To avoid a possible ambiguity in the latter meaning we shall write Ci{n). For 
k = only one root exists and is equal to unity. In the limit k — > the other 
two roots explode to infinity as ±k . The two branches (ci and C2) merge for 
«2 = \/2/27. The branch C3 is negative for k € (0, 1) and for k > K2 it is the 
only root of the relation ([T4"f . It is irrelevant in what follows and will not bother 
us any more. If k > 1, the equation has no real roots. Note, that parameter 7 
does not enter into the relation (TH|) . 

A linear perturbation around any of the constant solutions obeys the differential 
equation 



5 The Q-balls 



g = 




+ \/l + n 2 g 2 . 



(14) 



g'{ + -g[ + m\gi = 0, 



(15) 
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Figure 3: The three roots of the relation (IT?)) . 



where 



, (1 + c 2 k 2 )(3 + c 2 K 2 + c K 2 VT 



m = 7—^ - — ; : -• (16) 

1 + c 4 k 4 + c 2 k(2k-7) v ; 

c stands for a value of the perturbed constant solution. Roughly speaking, the 
solutions of (fT5j) oscillate with decreasing amplitude if m 2 > 0, otherwise they 
explode. In the expression (|16p the parameter 7 plays an important role as it 
makes the denominator vanish for some values of k. The infinite mass corre- 
sponds to a constant solution of the equation © being simultaneously a root 
of the relation ([7]). In this case the equation (|15D is not the correct one. To 
examine, how solutions from its close neighborhood behave, an analysis of the 
equation © with a-i — is required. 

Figure U shows how m 2 changes with n. It makes clear, that for small 7 the 
infinite mass appears for the solutions from the branch c\. The bigger 7, the 
smaller value of c\(k) for which the linear approximation fails. Finally, for 
7 > 4\/6/5 the infinite mass appears on the lower branch C2. When increasing 
7, the denominator in (|16|) vanishes for smaller and smaller n. 
Now we can describe the Q-balls present in the model. To this end we concen- 
trate on a few generic profile functions. Let us first consider the case of small 7. 
To be more specific, we assume that C2(n) <C <7r( K o)- It means, the structure 
described in the third section does not influence the solutions. We build the 
approximate Q-ball profile functions as a linearly perturbed C2(k) solution. It is 
possible if m 2 > 0. In analogy to the signum-Gordon model the solution reads 



9(y)=< \ Mvo)J , (17) 

[0 my>y a 

where the point yo w 3.8317 is the smallest y > satisfying the condition 
Jq(Ho) = and m — V to 2 . The fact, that the term m 2 weakly changes with /?, 
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Figure 4: Exemplary graphs showing the dependence of m 2 defined by for- 
mula (|16l) on k for a given 7.7 = 0.5 in the left panel, 7 = 3.5 in the right one. 
to 2 = marks the merging point of the branches c\ and c 2 . 



manifests itself in a scaling of the width of the solutions, since 

?/o = mur Q ur , 

where tq is a physical, not rescaled variable. Hence, tq ~ uj^ 1 . Such scaling was 
reported for spinning Q-balls in pQ. 

The approximation has an advantage: it gives hints on the regions in the pa- 
rameter space, where no Q-ball profile function is possible. The approximation 
fails if the starting point g(0) — c 2 (l — Jo(0)/Jo(j/o)) or "ending point" g = 
do not belong to the attraction region of the perturbed solution. Obviously it 
is the case, if g(0) > c\(k). This condition makes the range of allowed k values 
shrink to (0,0.186), where 7 is set to zero. 

For k close to zero the described approximation works well. The constant solu- 
tions c\ and C3 are too far to anyhow modify the perturbation. This is why the 
function (fT7|) is fairly close to the numerical solution, see figure El The bigger k 
the worse the approximation works. The explanation is rather easy: for bigger k 
the solutions start from a point in a close neighborhood of the constant solution 
c\ with to 2 < 0. The profile functions stay then almost unchanged on a rather 
long interval, see figure [5] Eventually, having left the vicinity of C\ they oscillate 
around the solution c 2 . Matching with the vacuum value is possible if K is not 
too big. Taking too large value of k gives rise to a function oscillating around 
c 2 with an amplitude too small to reach the vacuum value. 
Let us now move to larger values of the parameter 7. For small k both the 
constant solution Ci(/c) and the value of g r are again irrelevant to a small per- 
turbation around C2(k). Thus, the approximation (IT71) is fairly close to the 
actual solutions. For larger k the Q-ball profile functions start at the origin 
with a value close to a root of the relation ([7]). This results in sharpening the 
maximum at the origin as described previously. Thus, we obtain peakonlikc 
solutions, see figure [SJ It seems, that a maximal K admitting a Q-ball profile 
function emerges as a result of mutual approaching c\ and c 2 with increase of 
k, analogously to the above described case of large k and small 7. 
We assume, that in our model there is a region of small values of the parameter 
k for each value of 7, where the approximation (|17[) works well. This is backed 
by the figures 1 and 3. They show, that for n small enough both the constant 
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Figure 5: The examples of the Q-ball profile functions. In the upper panel 
7 = 0.5, on the left picture k = 0.03, on the right n — 0.185. In the lower 
panel 7 = 55, k = 0.0001 on the left and k = 0.00083 on the right picture. The 
approximate solutions are drawn with green lines, the exact numerical solutions 
with the black ones. 



solutions ci and the roots of the equation ([7]) are far from C2 and they are 
unlikely to influence a perturbation around this solution. This region shrinks 
when increasing 7 but probably never vanishes. It is noteworthy, that in case of 
a regular potential there is a maximal oj (minimal k) allowing for the appear- 
ance of Q-balls. This is why we expect the solutions not to exist for arbitrary 
large values of 7 then. 

There has been some interest in the so-called extreme baby Skyrme model, see 
e.g. [5], [TO]. It is a theory defined by a modification of the Lagrangian (TTJ). The 
modification consists in omission of the first term in ([T]). This theory is seen as 
a limit /3 — > 00 in the baby Skyrme model. For a particular potential it has been 
proven that there are no Q-balls in the extreme baby Skyrme model, see [lOj . 
The just described mechanism of expulsion of the Q-balls from the theory with 
increase of 7 suggests that there are no Q-balls in the extreme Skyrme model, 
no matter what potential is taken into consideration. 

6 Conclusions 

In this paper we have considered in detail how the additional terms modifying 
the Laplace operator present in the baby Skyrme model influence the equation 
for the Q-ball profile function. These terms are responsible for appearance of 
peakons in the investigated theory. As we have mentioned above, it would be 
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interesting to study the impact of this modification in a broader perspective. 
Next we have demonstrated the failure of the power expansion in terms of a small 
parameter in the theory with nonsmooth field force. The parameter dependence 
of the matching point may be blamed for this. Another way to characterize the 
parameter dependence has been devised. In this schema the Q-ball solutions 
appear as relevant perturbations of a constant solution of the equation for the 
profile function. Such a simplification of a fairly complicated equation bears not 
only qualitative understanding of the domain of existence of Q-balls, but also a 
good quantitative description of the solutions in some range of parameters. 
The work may be continued in many directions. Most obvious is a detailed 
investigation of the spinning Q-balls and gaining some qualitative (and quanti- 
tative if possible) understanding of these solutions. The role Q-balls may play 
in the baby Skyrme model is also not known. To this end their stability should 
be studied. It would also be interesting to check if similar problems appear in 
the topologically nontrivial sectors. 
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